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We classify the local bifurcations of one dov quantum billiards, showing 
that only saddle-center bifurcations can occur. We analyze the resulting 
planar system when there is no coupling in the superposition state. In so 
doing, we also consider the global bifurcation structure. Using a double-well 
potential as a representative example, we demonstrate how to locate bifurca- 
tions in parameter space. We also discuss how to approximate the cuspidal 



loop using AUTO as well as how to cross it via continuation by detuning 
the dynamical system. Moreover, we show that when there is coupling, the 
resulting five-dimensional system — though chaotic — has a similar underlying 
structure. We verify numerically that both homoclinic orbits and cusps oc- 
cur and provide an outline of an analytical argument for the existence of 
such homoclinic orbits. Small perturbations of the system reveal homoclinic 
tangles that typify chaotic behavior. 

MSC NOS 37N20, 37K55, 37M20 
1 Introduction 

Quantum chaos is an underdeveloped area of dynamical systems theory. One 
purpose of studying it is to generalize the notions of classical Hamiltonian 
chaos to the quantum regime, which has not yet been done in a universally 
accepted manner. One type of quantum chaos is often called semiquantum 
chaos, as these systems consist of classical (Hamiltonian) variables coupled 
with quantum variables. |Q Vibrating quantum billiards are a representative 



example of semiquantum chaotic systems. ||17| , pi[] They may be used as mod- 
els for quantum well microdevice components (such as quantum dots and 
quantum wires), Fermi accelerators]^, and intranuclear particle behavior. 

In the present paper, we consider the bifurcation structure of vibrating 
billiard systems. Quantum billiards describe the motion of a point particle 
undergoing perfectly elastic collisions in a bounded domain. The particle's 
motion is described by the Schrodinger equation with Dirichlet boundary 



conditions. One defines the "degree-of- vibration" (dov) of a billiard as the 
number of boundary dimensions that vary with time. If the boundary is 
time-independent, the billiard is said to have zero dov. The linear vibrating 
billiard and the radially vibrating spherical billiard have a single dov, and 
the rectanglular billiard with time-varying length and width has two dov. 

A zero dov quantum billiard exhibits only integrable behavior if it is 
globally separable. [|I^ It must, for example, be convex and describable by 



a non-composite geometry. ||15|, [T^ If part of the billiard is concave (or is 
composite, like the stadium billiard), it may behave chaotically, as it shares 
many of the instability properties of the Anosov diffeomorphism. Bliimel 
and EsserS found quantum chaos in the one-dimensional vibrating quantum 



billiard. Liboff and Porter |jT^ extended these results to spherical quantum 



billiards with vibrating surfaces and derived necessary conditions for chaotic 



behavior. They also generalized their results to other one dov billiards. |l21 
The purpose of the present paper is to examine bifurcations in single dov 
quantum billiards that occur when one alters the potential in which the 
billiard resides. 

Vibrating quantum billiards are important for several reasons. Though 
an idealized model, they are nevertheless useful for the study of quantum 
chaos. From a more practical standpoint, vibrating quantum billiards may 
be applied to several problems in physics. The radially vibrating spherical 
quantum billiard, for example, may be used as a model for particle behavior 
in the nucleusl^^ as well as for the quantum dot microdevice component. ]T5 . 



Additionally, the vibrating cylindrical billiard may be used as a model for 



the quantum wire, another microdevice component. [|25[ Other geometries of 



vibrating quantum billiards have similar applications. They may also be used 
as models of Fermi acceleration of cosmic rays. |^ The study of quantum chaos 
in vibrating billiard systems is thus important both because it expands the 
mathematical theory of dynamical systems and because it can be applied to 
problems in nuclear and mesoscopic physics. 

In the present paper, we show that saddle-centers are the only type of 
bifurcations that can occur in one dov quantum billiards. When there is 
no coupling in the superposition state, we show how to analyze the resulting 
planar system analytically and numerically. Considering a double-well poten- 
tial as a representative example, we demonstrate how to locate bifurcations 
in parameter space. We also discuss how to approximate the cuspidal loop 
using AUTO as well as how to continue past it by detuning the dynamical 
system. We also mention a shooting method for a more detailed analysis 
of the cuspidal loop. Moreover, we show that when there is coupling, the 
resulting five-dimensional system — though chaotic — has a similar underlying 
structure. We verify numerically that both homoclinic orbits and cusps oc- 
cur, and we outline an analytic argument that demonstrates the existence 
of homoclinic orbits. Small perturbations of the system reveal homoclinic 
tangles that typify chaotic behavior. 

2 Equations of Motion 

The goal of the present project is to examine the behavior of one dov quan- 
tum billiards in various potentials in order to determine the effects of the 
potential on the behavior of the system. In particular, we analyze bifurca- 



tions of equilibria both analytically and numerically. We consider a two-state 
superposition solution to the vibrating billiard, and we examine the above 
problem for both the case in which the two states experience coupling and 
that in which they do not. 

The present problem is described by the Schrodinger equation with solu- 
tions that are constrained to vanish on a time-dependent boundary a(t).|]2l[| 
That is, 

^(r, t; a{t)) = for r = a{t). (1) 

Because of the time-dependence of the boundary, the above boundary condi- 
tion is nonlinear. The (mathematical) problem at hand is to find a boundary 
a{t) such that Dirichlet boundary conditions are satisfied on it. One can 
then, in principle, insert a{t) into the eigenf unctions in our normal mode ex- 
pansion of the wave i^{r,t] a{t)) in order to obtain nonlinear normal modes. 
When the radius a{t) behaves chaotically, the nonlinear normal modes are 
examples of quantum-mechanical wave chaos. One derives coupled nonlinear 
ordinary differential equations for a{t) (and other variables) using a Galerkin 



approximation. |TT|, |22[ Considering a two-term superposition state then cor- 
responds to taking a two-term Galerkin projection. The equations of motion 
that one obtains depend on whether a particular subset of the quantum num- 
bers in the two states are the same.||2T| For the case of the radially vibrating 
sphere, the quantum numbers in question are the orbital and azimuthal quan- 



tum numbers / and m, respectively. [|T^ 

If these quantum numbers are the same in the two states, there is a 



coupling between them. The evolution of the system is then described by 
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and 



" = Ji' (^) 



oa a^ 



where x, y, and z are Bloch variables^ 



a; = Pl2 + P21 , 1/ = i{p2l - Pl2),Z = P22 - Pll , 

^2+ 1/2 + ^2^1^ (7) 



Pmn = AmA^ is the density matrix[|l6|, a is a displacement, P is its conjugate 



momentum, M is the mass of the billiard, m <^ M is the mass of the confined 
particle, /i > is the coupling coefficient between the two eigenstates, V = 
V{a) is the potential of the billiard boundary, uq = ^^^, e± = 2^^^ ^^'^ 
61 and 62 (62 > ei) are the energies of the two states. It has been shown that 
these equations exhibit quantum chaotic behavior, pij 



If there is no coupling between the two eigenstates, the evolution of the 
system is described by a one degree-of-freedom Hamiltonian. The equations 



of motion are 

P 

M' 

where 

\ = 2{e,\Ci\^ + ei\C2\'), (9) 

Ci and C2 are constants such that |Cip+ IC2P = 1. The energy parameter A 
is necessarily positive because ej > and the |Cip correspond to probabihties. 

3 Integrable Case: Absence of Coupling 

For the planar case, all (a, P) that satisfy a = P = are equilibrium points. 
Each one is of the form (a, 0), where a satisfies 

f^(«,o) = 4- (10) 

The eigenvalues of the integrable system (||, ^ at the stationary point (a, 0) 
are 



These eigenvalues are either real with opposite sign or are pure imaginary 
pairs, so in the linear analysis, each equilibrium is either a center or a saddle 
point. If 



then every equilibrium point is a linear center. This holds, in particular, when 
the potential V{a) has a single minimum (single- well potentials). Previous 
studies have focused on the harmonic potential 

V{a) = ^ia-aof = V2{a-aoY. (13) 

Another interesting single-well potential is the quartic potential 

V{a) = ^{a-ao)'^V,{a-ao)\ (14) 

% 

In the above equations, the Vj are dimensionless parameters. 

It is insightful to examine the evolution (particularly in the chaotic case) 
of the quantum billiard system with the above two potentials and derive 
a mechanical anology in terms of spring stiffness. Given the same initial 
conditions and the quartic and quadratic potentials above (and assuming 
V2 = V4 for ease for comparison), one observes that the phase-plane trajectory 
described by the evolution of a and P for the quartic potential (|1^) has 



a larger radius of curvature (that is, a smaller curvature). For all initial 
conditions, the trajectory in the quadratic potential has a larger maximum 
a. For initial conditions with sufficiently small a(0), the quadratic potential 
induces trajectories with smaller maximum |P|, but the quartic potential 
eventually gives a larger maximum |P| as a(0) is increased. 

An equilibrium for which A < is a saddle point. Since the present 
system is a Hamiltonian system, it is invariant under reflection about the 
a-axis, so that the eigenvectors representing the local stable and unstable 
manifolds are mapped to each other under this reflection. Since the only 
other possible types of equilibria are centers, it follows that if there is at 
least one saddle point, the system must have saddle connections. If there is 
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exactly one, the connection is a homoclinic orbit, and there must be two of 
them emanating from the saddle because the system must have a center (in 
the right-half plane) on each side of the saddle. 

As one considers increasingly excited states of the system (corresponding 
to larger quantum numbers), the energy parameter A is increased. Each sad- 
dle will eventually become a center as this occurs. The quantity A vanishes 
at such an equilibrium point. The stationary point then has a double zero 
eigenvalue with the Jacobian 



(15) 



so that the saddle-center bifurcation that occurs has codimension two and 




gives rise to a global bifurcation corresponding to the breaking of the separatrix. |]I4 



111 To find the conditions satisfied at this point one can either solve the simul- 
taneous equations a = 0, cr = or find the points at which the Hamiltonian 
has a double zero (which is equivalent to solving the system of equations 
H{d, 0) = 0, ^ (a, 0) = 0). If the potential has a constant term Vq, it does 
not change the evolution of the system since it does not appear in the equa- 
tions of motion. We thus let Vq = without loss of generality. One finds 
that a saddle-center bifurcation occurs when 

r)V 
A = #^(a,0) (16) 

at the point (a, 0) satisfying 

dV , , ad'^V , , 

Any solution to (|l^ giving A < is discarded as nonsensical. 
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At the saddle-center bifurcation point, the stable and unstable eigenvec- 
tors of the equilibrium coincide along the a-axis, so that the stable and unsta- 
ble manifolds overlap near the stationary point. This cusp causes difficulties 
in numerical continuation attempts, as standard continuation techniques fail 
for this bifurcation. One observes that the two homoclinic orbits that exist 
when A < have coalesced into one. (As A increases, the homoclinic orbit 
on the left shrinks, becoming a single point at the saddle-center. The orbit 
has infinite derivative with respect to arclength at the saddle-center point.) 




Figure 1: Approximate cuspidal homoclinic orbit. The initial point used was 
(0.7886751,0.001), just above the equilibrium. 



As a specific example of this phenomenon, we consider the quartic poten- 



tial 



V{a) = Vi{a - ao)^ + V^ia - aof + ^2(0 - oq)^ + Vi{a - ao), (18) 

where Vi = ^. This potential has either one well or two. In the latter 
case, there is a single saddle-center bifurcation point. One can find a and 
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A exactly in this case, since the former is a root for a degree-three polyno- 
miaL We initially restrict ourselves to the case in which V3 = 14 = 0, since 
all the dynamics of interest remain in this simpler case. Note that the po- 
tential is symmetric about oq. For numerical purposes, consider the special 
case corresponding to the parameter values ao = 1,V2 = —1, and V4 = 1. 
There is a saddle-center bifurcation at A = 0^2 [3 -|- a/s] v^ ~ 0.1888176. 
The corresponding stationary point is (| + g jO) ~ (0.7886751,0). Using 
DsToolfT^, we plotted an approximation of the homoclinic orbit emanating 
from this equilibrium (Fig. |l|). 

More generally, one can consider even polynomials of higher degree in 
order to examine vibrating quantum billiards in an A^-well potential. If the 
polynomial is of degree six or higher, one may not be able to solve for a exactly 
in terms of radicals by Galois theory 0, since equation (|1^ is polynomial of 
degree at least five. If its degree is exactly five, one can solve for a exactly 
in terms of elliptic functions. [^ 



Consider the problem of starting at A < A and attempting to continue 
along the bifurcation curve past the saddle-center. Using AUTO[0, H and 
the homotopy method encoded in HomCont, we followed the two homoclinic 
orbits for A = 0.15 (Fig. ^. The saddle connection in the present system 
has a codimension greater than one, as both regularity and non-degeneracy 
conditions are both violated. 0, |^ The present situation is degenerate because 
for all A < A, there are two homoclinic orbits emanating from the saddle 
point. Regularity is violated because the saddle point's two eigenvalues are 
negatives of each other. (Moreover, their eigenvectors are related by reflection 
across the a-axis, since the system is Hamiltonian.) 
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Figure 2: Homoclinic orbits emanating from (0.8916637, 0) for A = 0.15. The 
label 12 refers to the right homoclinic orbit, and the label 6 refers to the left 



one. 



Because the present system is degenerate and irregular, one cannot con- 
tinue (in A) past the saddle-center directly, as described in the AUTO manual. [ 
Hamiltonian systems possess a continua of homoclinic orbits, and there are 
numerical schemes that allow one to handle this phenomenon. One can ex- 
ploit the reversibility of the system by computing only half of a given saddle 
connection. However, the cusp at the saddle-center point prevents this from 
working for the present system. AUTO is incapable of continuing past a 
cusp, because ^ vanishes there and the pseudo-arclength step must become 
vanishingly small for such a continuation step to be successful. Because of 
machine precision, this cannot occur, and so one must "detune" the system 
to continue across the cusp. 

In general, Hamiltonian systems are described by 



X 



JVH{x,X), X e 



h2n 



(19) 
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where 



J^l ° ^1 (20) 

-/ 



is the canonical 2n x 2n symplectic matrix |^. One can detune the sys- 
tem by using a small perturbation parameter e as follows. ^ The perturbed 
dynamical system, 

x = JVH{x,X) + eVH{x,X), (21) 

is no longer Hamiltonian, but the perturbation was constructed so that the 
locations of all equilibrium points are preserved. With this detuning, the 
saddle-center bifurcation becomes a saddle-node bifurcation (the eigenvalues 
of the stationary point are now of the form a ± y/F (a j^ 0) rather than of the 
form ±v&). One can then continue A past this point using AUTO (without 
utilizing the HomCont package). Using this technique, one can compute the 
value of A at which the saddle-center bifurcation occurs as well as the cusp 
point of the homoclinic orbit corresponding to that value. Moreover, once 
one has successfully continued past the cusp, one can simply let e — > and 
thereby work with the system when A > A. This method of continuation is 
useful as long as one needs to get past the cusp rather than do computations 
at the cusp itself. For the present system, we used e < 0, since in that 
case the equilibria that are centers for e = become stable spirals. The 
continuation curve (in A) is shown in Figure ^. 

In general, AUTO has difficulties near cusps. As with DSTool, one can 
approximate the cuspidal homoclinic orbit using AUTO. In order to do this, 
one provides initial values for the HomCont continuation parameters (corre- 
sponding to the initial point in the (a, P)-plane) to the right of the saddle 

13 



0.70 0.90 

0.80 larrbda 



Figure 3: Continuation of the detuned system in the parameter A. 

point (a*,0). This allows AUTO to continue along the homoclinic orbit for 
values of A closer to A than if one had started as close to the saddle point as 
machine precision would allow. For the present example, the closest accurate 
plot we obtained was for A = 0.1887. The right homoclinic orbit is pictured 
in Figure § and the left one is pictured in Figure §. Observe that the one 
on the right looks like it has a cusp at the saddle point and that the left one 
is very small. As the saddle-center is approached, the left homoclinic orbit 
shrinks to a single point and the right one becomes more cusplike. 



4 Chaotic Case: Presence of Coupling 

There are other methods one can use to perform analysis near the cusp. 
One can, for example, use shooting methods. pO|] (AUTO uses a collocation 
method.) Saddle-centers are a degenerate case of Takens-Bogdanov (TB) 



bifurcations ||Tl[, so one can add dummy parameters and analyze the cuspidal 
loop by computing the locus of a TB bifurcation in parameter space, moving 
along the TB curve until one finds a Hopf bifurcation of another equilibrium. 
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Figure 4: Right homoclinic orbit for A = 0.1887. 

and then following the evolution of the periodic orbit created in the Hopf 
bifurcation as the parameters follow the TB bifurcation curve. If a cuspidal 
loop exists, this method will find it when the periodic orbit collapses into the 
cusp point. 

For the chaotic case, the equilibrium points satisfy x = P = y = 0, 
z = ±1, and ^ = ^(e+ ± e_), where the factor of ±1 in the last quantity 
corresponds to the equilibrium value of z. The eigenvalues for the equilibria 
of the chaotic system (Eqs. ^|-^ are of the form 



a = 0, ± 



-2M{a ± y/P) 



(22) 



2Ma2 

where both a and (3 have terms whose signs depend on whether z is 1 or —1. 
The quantities a± and f3± are given by 



a+ = a^^^v + u^M - Ve_ + 6(e+ + e_), 



(23) 
(24) 



13, 



I, + I+ + I+ + I+, 



(25) 
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where 



Ii = a'' 



da^ 



It = -K-Y (I2a^e+ - 2ujlMa^ + 12a\_ - 8^^a^e_) , 
J3+ = 16/i^e^ - Scj^M/i^ _ 48/1^6^ - 12cj^Me+ + 72e+e_ - 12a;^Me_, 

/4+ = -48fih+e- + cu^M^ + 36(e^ + e^), (26) 

and 



P.=h + I,+I,+I,, (27) 



where /i is as before, 



I^ = -— (I2a^e+ - 2a;2Mo^ - 12a^e_ + 8^^a^e_) , 

/g- = 16fi^e^_ + Su^Mfi^ - 48^2e! - 12a;^Me+ - 72e+e_ + 12cjgMe_, 

/4~ = 48^2e+e_ + co^M^ + 36(e^ + e^). (28) 

Analogous to the planar system, only a generalization of saddle-center 
bifurcations can occur. As the energy is increased, a bifurcation corresponds 
to an increase in the dimension of the center manifold by two as a pair of real 
eigenvalues of opposite signs becomes a pair of pure imaginary eigenvalues. 
For a double-well potential, the only possibility is a jump in the dimension of 
the center manifold from three to five. For more complicated potentials, there 
may be parameter values with a one-dimensional center manifold. As before, 
one can determine the location of this bifurcation by finding the equilibria 
for which H{a, P) has a double root. One again finds that the equilibrium 
point (0, 0, ±1, a, 0) at the bifurcation satisfies 

dV ad'^V 

16 




0.775 0.780 0.785 0.790 0.795 0.800 0.805 



Figure 5: Left homoclinic orbit for A = 0.1887. 

Finding the parameter values at which this occurs is not the simple issue 
it was in the previous case. In the planar case, A was a probabilistic weighting 
of two energies, so it could be varied continuously past the bifurcating value 
A. However, the parameters ei, £2 in the present system are part of a discrete 
energy spectrum, and so one cannot vary them continuously. In practice, 
therefore, this model predicts superposition states on each side of the bifur- 
cation, but one does not expect to observe the system at a bifurcating value 
of (ci, €2) because the set of all (ei, 62) that correspond to bifurcating values 
has measure zero. (One could vary V2 continuously if one wanted to examine 
bifurcations corresponding to a change in the quartic potential. If V2 is nega- 
tive and sufficiently small for a given V4 or if it's positive, the system will not 
exhibit a saddle-center.) Additionally, numerical observations indicate that 
the bifurcation occurs at low energies (corresponding to superpsoitions of low 
quantum number states), so that for a given billiard system, most superposi- 
tion states will exhibit an evolution with a five-dimensional center manifold. 
The bifurcation under study may thus occur as one considers superpositions 
of increasingly excited states of the quantum billiard. 
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Numerical investigations have shown that homochnic orbits must exist 
for this five-dimensional system (Fig. |^). Slight perturbations away from the 
homoclinic orbits lead to homoclinic tangles (Fig. |^, which are traditionally 
analyzed using symbolic dynamics. The details have not been worked out, 
but the existence of a homoclinic orbit for the present case can be shown an- 
alytically as follows. There is only one saddle point, so any saddle connection 
would have to be a homoclinic connection. By Hamiltonian symmetry and 
the existence of a center to the right of the "saddle-like" (in the sense that it 
has one-dimensional stable and unstable manifolds) equilibrium, there must 
exist a structure to its right that looks like a homoclinic orbit when projected 
into the (a, P)-plane. It may not be a homoclinic orbit, because one must 
consider the value of [x, y, z) at the point (a, 0) where the projection inter- 
sects the a-axis. One thereby considers the 2-sphere S"^ and two trajectories 
on it that start at the same point. One curve begins at t = oo from the stable 
manifold, and the other starts at t = — oo from the unstable manifold. One 
looks at the intersection of one trajectory at time t with the other at time 
—t. If one can prove that such a point exists at a time T, then repeating the 
argument shows that such an intersection occurs at infinitely many times. 
To complete the proof, one must show that P = at one of these points. 

Note that the homoclinic orbit for the chaotic system (Eqs. 0-^) em- 
anates and terminates from a nonhyperbolic equilibrium, which increases 
the difficulty of numerical studies, p, ^ As with the planar system discussed 
earlier, there is also a cusplike structure as the stable and unstable manifolds 
coalesce along the a-axis. 

As an anology, consider two undamped springs, one with a linear restoring 






a 








Figure 6: Poincare section projected into the (a, P)-plane demonstrating that 
there must exist a homochnic orbit. 

force {Fi = —kx) and a stiffer spring with a cubic one (F3 = —kx^). These 
two springs (with mass m = 1) are described, respectively, by the differential 
equations 

x + kx = 0, (30) 

and 

x + kx^ = 0. (31) 

One observes the same stiffness phenomena for these spring systems as we 
did for the integrable case of the vibrating billiard. (For example, the phase 
space trajectory of the cubic spring (Eq. ^) has a larger radius of curvative 
than the analogous one in the linear spring system. There is a correspondence 
with the other properties we discussed as well.) By analogy with mass-spring 
systems, we thus conclude that it is reasonable to consider the "stiffness" of 
the potential in which a quantum billiard resides as an object of interest. 
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Figure 7: A homoclinic tangle projected onto the (x, y)-p\a.ne of the Bloch 
sphere. 

The analogy with spring stiffness carries through in the chaotic case as 
well as in the planar case. For the same initial conditions, we consider the 
potentials 



V{a) = Viia - aof 



(32) 



and 



V{a) = V,{a-ao)\ 



(33) 



where for ease of comparison, V2 = V4. As shown in Fig. p, the quartic 
potential gives trajectories with a larger radius of curvature than those in 
the quadratic potential. It thus makes sense to consider a potential's stiffness 
in the chaotic case as well as in the integrable one. 
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Figure 8: Chaotic Poincare maps in the (a, P)-plane for bilhards in both 
a quartic and a quadratic potential The plot from the quartic potential 
is the one with a smaller maximum value for the distance a. Notice that 
the trajectory corresponding to the quartic potential generically has a larger 
radius of curvature. 

5 Conclusions 

We showed that saddle-centers are the only type of bifurcations that can 
occur in one dov quantum billiards. When there is no coupling in the su- 
perposition state, we showed how to analyze the resulting planar system 
analytically and numerically. Considering a double-well potential as a rep- 
resentative example, we demonstrated how to locate bifurcations in param- 
eter space. We also discussed how to approximate the cuspidal loop using 
AUTO as well as how to continue past it by detuning the dynamical sys- 
tem. We also mentioned a shooting method for a more detailed analysis of 
the cuspidal loop. Moreover, we showed that when there is couphng, the 
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resulting five-dimensional system — though chaotic — has a similar underlying 
structure. We verified numerically that both homoclinic orbits and cusps oc- 
cur. Small perturbations of the system reveal homoclinic tangles that typify 
chaotic behavior. 
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Figure Captions 

Figure 1: Approximate cuspidal homoclinic orbit. The initial point used was 
(0.7886751,0.001), just above the equilibrium. 

Figure 2: Homoclinic orbits emanating from (0.8916637,0) for A = 0.15. 
The label 12 refers to the right homoclinic orbit, and the label 6 refers to the 
left one. 
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Figure 3: Continuation of the detuned system in the parameter A. 

Figure 4: Right homochnic orbit for A = 0.1887. 

Figure 5: Left homochnic orbit for A = 0.1887. 

Figure 6: Poincare section projected into the (a, P)-plane demonstrating 
that there must exist a homochnic orbit. 

Figure 7: A homochnic tangle projected onto the {x, ?/)-plane of the Bloch 
sphere. 

Figure 8: Chaotic Poincare maps in the (a, P)-plane for biUiards in both 
a quartic and a quadratic potential The plot from the quartic potential 
is the one with a smaller maximum value for the distance a. Notice that 
the trajectory corresponding to the quartic potential generically has a larger 
radius of curvature. 
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